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The paper presents results on the stability and dynamic behavior 
under large signal conditions of nehvorks consisting of transistors and 
sources connected to a linear, passive, memoryless subnetwork. The tran- 
sistors' model incorporates various nonlinearities. A characteristic common 
to the main results of the paper is that they relate to properties of the tran- 
sistors alone and, hence, are independent of the passive part of the network. 

Sufficient conditions are obtained for asymptotic and bounded input- 
bounded output stability. The conditions impose restrictions on some of 
the physical constants of the transistors' model. These conditions have 
an interesting physical interpretation in terms of temperature differentials 
in the transistor junctions. In particular, any transistor with the expo- 
nential type of static diode characteristic is passive only if the ratio of the 
junction temperatures lies inside an interval determined by the a's. 

In the state space of the network there exists a well-defined region R 
specified by the transistors' model with the property that constant terminal 
stales in R are independent of initial conditions. The region R is in a 
certain sense maximal. 

I. INTRODUCTION AND DERIVATION OF THE DIFFERENTIAL EQUATION 

The network considered is shown in Fig. 1 and it consists of a number 
of transistors connected to a subnetwork composed of voltage sources, 
current sources, and linear, passive, memoryless devices such as resistors 
and transformers. Sandberg 1 has analyzed the dynamic behavior of such 
networks and has denned a class which exhibit various features of 
stability. The results of this paper are essentially different in that they 
relate to properties of the transistors alone and, hence, are independent of 
the passive part of the network. 

The analysis is concerned with certain network-theoretic properties of 



* This paper was presented at the 1970 IEEE International Symposium on Circuit 
Theory, Atlanta, December 1970. 

f The sequence of names was chosen by flipping a coin. 
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Fig. 1 — General network containing transistors, sources, resistors, and transformers. 



transistor models that are assumed in extensively used Network-Analysis 
Programs, such as NET 1, CIRPAK, CIRCUS, etc., (see Ref. 2). The 
intent of the paper is to study the large signal behavior of transistor 
networks using this model and certain generalizations of it. Each transis- 
tor is represented by a model of the type shown in Fig. 2 which takes into 
account the nonlinear dc properties as well as the presence of nonlinear 
junction capacitances. Six parameters are associated with the model: 
a,, a 2 , T], t 2 , Cj, and c 2 , all of which are positive; also a,, a 2 < 1. The two 
nonlinear static diode functions are denoted by /,(•) and / 2 (-)- Initially 
it is assumed that f x (-) and / 2 (-) are monotone, strictly increasing 
mappings of the interval (— °o, °o) into itself; /i(0) = / 2 (0) = 0, and 
/,(•) and / 2 (-) are continuously differentiable on (— go, co). Further 
assumptions are made about fi(-) and / 2 (-) in the course of the paper. 
Suppose the network has n transistors and with the polarity indicated 
in Fig. 2 let v 2i -i and v 2i respectively denote the emitter-to-base voltage 
and collector-to-base voltage of the ith transistor. Similarly, let * a «-i 



o 



a,U r 
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-► If = f , ( V| ) 
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a 2 lf 



Lr = f 2 (V2)- 



A/vV 




Fig. 2 — Transistor model. 
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and i 2 i respectively denote the emitter current and the collector current 
of the ith transistor. An identical scheme of subscripting is used to 
define /,-(•)> c,, and t, for i = 1,2, ■ ■ ■ , 2n. v, i, and F(-) are 2»-column 
vectors formed by arranging {»,•}, {*«}, and {/,•(•)} respectively. In 
applying the current law to the transistors' model it follows that 

i = f t C(v) + TF(v) (1) 



where 



[C(v)] t = c iVi + r f /,(t,) 



and T = T l T 2 • ■ ■ © T n the direct sum of n 2 by 2 matrices 7\ in 
which 



T, = 



i -«; 

L-«i i . 

for .7 = 1, • • • , n. 

The subnetwork to which the transistors are connected is assumed to 
impose a constraint of the form 

i = -Gv + b (2) 

in which G is the conductance matrix and hence G ^ 0,* and b is an 
element of all real bounded continuous 2;?-vector-valued functions of t 
on [0, «>). 

From (1) and (2), 

| C(v) + TF(v) + to = b. (3) 

The above equation is also derived in Ref. 1. Since all of the c, and t, 
are positive and each of the /, ( • ) is continuous and monotone increasing, 
C~'(-) exists; obviously C~'() is also strictly monotone increasing 
and C"'(0) = 0. 

II. ASYMPTOTIC STABILITY FOR THE UNFORCED SYSTEM 

In this section we show how to prove an intuitively reasonable result 
concerning the asymptotic stability of the system. A Lyapunov function 
with a simple energy interpretation is introduced and used to prove 



*(?|0 indicates that G is a positive semidefinite matrix. Unless otherwise stated, 
G will only be assumed to be positive semidefinite. 
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stability. Let 

Uv) A Z f°' '' C7\a) do. (4) 

Since C,(-) defines the charge-voltage relationship in the nonlinear 
junction capacitances, L(v) is the total stored electrical energy in the 
network. It is easily verified that (i) L(v) ^ 0, equality holding only 
when v = 0, and (ii) L(v) — > °o as ||i»|| — > °o. Thus L(v) is a Lyapunov 
function if (d/dt)L(v) < 0. For the unforced network 

j L(v) = -v'TF(v) - v'Gv. (5) 

Clearly if v'TF(v) ^ with equality holding only if v = then L{v) is a 
Lyapunov function since G has been assumed to be positive semidefinite. 
It then follows from a well-known result in stability theory 3 that 
v(t) — » as t -* oo . 

Note that for constant y the term v l TF(v) expresses the net power 
flow into the transistors and if all the transistors are passive then 
certainly v'TF(v) > 0. However it is shown later that in certain abnormal 
conditions the transistors are not all passive and the judicious choice of 
the passive network allows one to extract energy from the transistors. 
To show this we first prove the following: 

Lemma 1: Given two real 2n-veclors x and y such that x'y > 0, there exists a 
G > for which the following holds 

y = Gx. (6) 

Proof. Define a 2n by 2n — 1 matrix Z such that the columns of Z span 
the (2n — l)-dimensional subspace orthogonal to x. Then Z l x = 0. 
Since by assumption y is not orthogonal to x, y is not an element of the 
range space of Z and, hence, the columns of Z together with y span E 2n . 
For these reasons the following construction of G, 

G-4-ftf + zz; (7) 

•*■ y 

suffices to prove the lemma. 

When for some Vy , v[TF(vi) < then by Lemma 1 there exists a 
G > such that 

TF(v t ) + Gv t = 0. (8) 

Hence with initial condition v(0) = v, the solution of the network 
equation (3) with b = is, 
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V(t) = V, 
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(9) 



and the system is not globally asymptotically stable. 

The above results may be summarized in the following: 
Theorem 1: If v'TF(v) >0V^0 then the network is globally asymp- 
totically stable. Furthermore, if there exists a v, such that v^Fipx) < 
then there exists a G > for xchich the system is not globally asymptotically 
stable. 

2.1 Sufficient Conditions for Asymptotic Stability 

Consideration is given to the positivity of v l TF{v). Due to the quasi- 
diagonal structure of T we need only consider in detail the behavior of 



<p(v 2i -i ,v 2i ) = (y 2 ,-i , p«) 



1 -a\ 
-at 1 






(10) 



For notational convenience a[, a\, v 2i -i, v 2i , / 2 ,-i(-)> and / 2 ,(-) are 
respectively denoted by a ly a 2 , v u v 2 , /i(-). and f 2 (-). 

Note that when v,v 2 g 0, <p(v) ^ and <p(v) = only if v, = v 2 = 0. 
Let Vi ?* 0, v 2 = pv u and p ^ 0. Then 

-fo) = vJM - ai»i/ 2 (t» a ) + v 2 f 2 (v 2 ) - ajWa/iCtfi) (11) 



and 



</(p) * 



= (l-p« 2 ) + / f ( gf(p-a I ). 



(12) 



If P = <*! then <?(p) > since 1 - «>a 2 > 0. If p ^ «i, fl(p) = if and 
only if 

1 — pa -i = fi{pV\) Q3) 

ai - P " /i(»i) 
(1 - pa 2 )/(ai - p) is plotted as a function of p in Fig. 3. The salient 
features of the function for p ^ are: 



1 ~ p " 2 ^ 1 for ^ p < a, 
a! — p <*i 



< for a, < p ^ — f • 
a 2 



< a 2 for — < p < °o 



(14) 



Now [/ a (p»i)//i(fi)] is a monotone, strictly increasing function of p and 



2840 



THE BELL SYSTEM TECHNICAL JOURNAL, OCTOBER 1971 
1-pa 2 




Fig. 3— Plot of (1 - pa 2 )/(<*i - p). 



[/2(pWi)//i(«i)]| P -o = 0. Clearly if 

Ufa) = «i 



and 



\« 2 



^ a 2 



V o, ^ 



(15) 



then no solution of (13) exists. Since g(p) is continuous in p, g(0) > 0, 
and g(p) ^ V p > 0, it follows that p(p) > V p ^ 0. Hence sufficient 
conditions for the positivity of v'TF(v) are: 



and 



Remark: 



/z,-i(y) ~ a! 



/«(») 



. V v ^ and z = 1, 2, • ■ • , n. 



Ui-i(alv) 



rr ^ «2 



(16) 



7 7T < for < p ^ a, 

> for — ^ p < oo 

<*i — p a 2 



■ V w ^ 



(17) 



is both necessary and sufficient for v'TF(v) > 0. However, the practical 
value of such a condition is limited. 
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2.2 Necessary and Sufficient Conditions for the Exponential Type of 
Nonlinearities 

The functions /,(•) in the conventional Ebers-Moll model and in the 
charge-control model 4 are of the following form: 

f 4 (v) = a, (exp X,y - 1). (18) 

Such exponential nonlinearities are subsumed in a class of nonlinearities 
with the following properties: 

(i) f.(y .) = a,/(X,?\), i = 1,2, ■■ ■ , 2n, where a, and X, are positive 

constants. 
(ii) /(•) is a monotone, strictly increasing function. 
(in) /(0) = 0. 
(iv) lim^-oo /(y) = — 1. 
(v) lim,-^ \j(pu)/f(v)] = for all p where < p < 1. 

When the transistors' static diode characteristics fall under this class 
then the following is true: 

necessary and sufficient conditions for vTF{v) ^ (v'TF(v) = 
only when v = 0) are 

a\ ^ a ^ - , ^ ^ A (19) 

a 2 i A;., a 2 

for i = 1, 2, • • • , n. 

Sufficiency follows from the definitions and (16). For the necessity 
part consider the case when v 2i -i — * — « with the remaining 2n — 1 
variables fixed. It is necessary that the coefficient of v 2i -\ in the expansion 
of vTF(v) approach a nonpositive value, i.e., 



lim \j 2i -M -a\f 2 ,(v)\ ^ (20) 



for i = 1,2, ■ • • , n. 
Similarly, 



lim |/„(») -ciU-M) S (21) 



for i = 1,2, • ■ ■ , n. Hence 

— a 2l -i + a\a 2 i ^ 



and, 



— a 2 , + a 2 a 2 ,_i ^ 
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or 

. ^ <hi-i ^ 1 



<*i ^ 



a 2i a 2 



(22) 



for t — 1,2, • ■ • , n. To obtain the remaining conditions consider equation 
(12) in the notation of this section. By property (v), 

lim |^M- = 
for all p such that p < -- 1 and i — 1,2, • • • , n. 

A 2 , 

Hence, a necessary condition for the positivity of v l TF(v) is: 1 — pa* 2 > 
is implied by p < X 2l _i/X 2l ; i.e., 

¥=* ^ A (23) 

A2i «2 

for i = 1, 2, • • • , n. Repeating the argument with the roles of v 2l _i and 
v 2i interchanged yields 

at g ^ (24) 

A 2 < 

for i = 1, 2, • • • , ». This concludes the proof. 

Remark 1: It is clear from the derivation that the conditions stated 
in equation (19) with the ratios replaced by the appropriate limits 
are necessary for the positivity of v'TF(v) when {/,(•)} are general 
monotone, strictly increasing functions. 

Remark 2: If condition (19) is violated then there exists b.v x , || t>i || < «> , 
such that v[TF(vi) < 0. The proof is as indicated below. 

Consider equation (12) and say Xi/X 2 > l/a 2 ; then for p such that 
l/a 2 < p < Xi/X 2 it follows from property (v) that there exists a V such 
that for Vi > V, g(p) < 0. Likewise, for a x /a 2 > l/a 2 and p(p > 0) such 
that (1 — «i(a 2 /ai)) + p{a 2 /a x — a 2 ) < ft-it follows from property (iv) 
that there exists a V such that for v x < V, g(p) < 0. For the remaining 
possibilities the proof follows by interchanging the roles of v x and v 2 in 
the definition of g(p) and proceeding as above. 

III. PASSIVITY 

This section is devoted to a discussion on condition (19) which has 
been shown to be both necessary and sufficient for the passivity of 



KT, and X " 


9 


A21-I 1-2 

x 2 , r, 
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transistors for which the static diode characteristics are of the form 

/,(»,) = o,(exp \iV t - 1) (25) 

in which a, and X, are positive constants. When the temperature in the 
transistor is uniform, then from the well-known physical model, 5 X, = 
(q/KT) where q, K, and T are respectively the electron charge, Boltz- 
man's constant, and temperature. However, when a temperature dif- 
ferential exists in the transistor junctions, for instance in a p-n-p 
transistor the temperature in the neighborhoods of the p-n and n-p 
junctions are respectively 7\ and T 2 , then the same physical model 
holds 6 and 

A.2.-1 = 

so that 

(26) 

Hence when the ratio of the temperatures at the base-collector and 
emitter-base junctions of at least one transistor, say the t'th, lies outside 
the interval [a\ , 1/aj] then there exists a 2n-vector v t such that v[ TF{v y ) 
< 0. Furthermore, by Lemma 1, there exists a positive definite matrix 
G such that the solution of the network equation (3) with 6 = and 
initial condition v x is v(t) m v x . Then, 

power delivered to thel = ^ = - v [ T F(v,) > 0. 
passive subnetwork J 

The temperature dependence of the o.'s may similarly be exploited to 
deliver power. The phenomenon of differential heating of transistors as 
described above is the basis of the thermocouple effect and thermo- 
electric generators. 5 

IV. BOUNDED INPUT — BOUNDED OUTPUT STABILITY 

It is shown in this section that for a positive definite matrix G, bounded 
inputs, and passive transistors there exists a bounded neighborhood of 
the origin which, loosely speaking, is sure to contain the forced response 
of the network. The norm used is the Euclidean norm. 

Now 

v l TF{v) > Vy^O (27) 
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and let 

II 6(0 II ^ K. (28) 

L(v) is as defined in equation (4). For the forced system, v ^ 0, 

j- L(v) = -v'TF(v) - v'Gv + v'b 

{XL 

^ t/~i I I * I. I 

< — V (rV + | V | 

^ —v'Gv -f- K || v || by Schwarz's inequality (29) 

^ -x„.i„ ii v ir + k ii v H 

where X min is the smallest eigenvalue of G and is positive by assumption. 
Hence when \\v\\ ^ K/X min , (d/dt)L(v) < 0. Thus there exists a T such 
that for t > T, 

IK*) II <r^- 

A m , n 
V. A PROPERTY OF TRANSISTOR NETWORKS 

Let R = {v'\ (v' - v)'T[F(v') - F(v)] > Vv ^ v'\. Then: (i) if G 
is positive definite, v' and b' are constant vectors such that TF{v') + 
GV = b', and »' e 72, then for all inputs which approach b' as t — * go , 
the state vector approaches y' independent of initial conditions; (ii) there 
exists a passive subnetwork for which the overall network is such that 
if the terminal state corresponding to an input approaching a constant 
vector is not in 7^, the closure of R, then there exists another input 
approaching the same constant vector for which the terminal state is 
different. 

Remark: It has been pointed out by Sandberg 1 that the independence of 
the steady state from initial conditions is a basic property, and in Ref. 1 
it has been proved that if (T, G) belongs to a certain class then the region 
R in statement i extends to the entire state space of the network. How- 
ever, the region R defined here is independent of G. 

In switching circuits with "memory" the dependence on initial 
conditions is a salient feature. Statement (i) observes that the design of 
such circuits be such that the steady states lie outside the region R, i.e., 
the steady state bias voltages violate the conditions which define the 
region R. Of course, the fact that the transistors are passive merely 
implies that R is non-empty, and their use in switching circuits with 
"memory" is not precluded. Statement (ii) states that in a certain sense 
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the region R is maximal. 
Proof: (0 

TF(v') + Gv' = b' (30) 

and 

| C(v) + TF(v) + Gv - b' =- dl \C{v) - C(v ')} 

+ T\F(v) - F(v')\ 4- <?(«-*') 

(31) 



(32a) 
(32b) 

(33) 



j C(v) + TF(T>) + Gv = 6(0 - V. (34) 







= Ht) - 


b' 


from (3) and (30). 


Define 








V 


= v - v' 






C(v) 


= C(v) - 


C(v'), 


and 










F(v) 


- F(v) - 


F(v'). 



From (31), 



Note that (?,(•) for i = 1,2, • • • , 2» are monotone, strictly_ increasing 
mappings of the real interval (— *>, oo) onto itself and C,(0) = 0; 
also there exists C7'(-)- Since v' t R, v'TF(v) > 0. Also by assumption 
for any e > there exists T such that 

|| 6(0 - 6' || ^ X min e for t ^ T (35) 

where \ rain is the smallest eigenvalue of G. On applying the results on 
bounded input-bounded output stability, as stated in Section IV, to the 
system described by equation (34) it follows that there exists a T' such 
that for / > T', 

II 3(0 || = 1 1 0(0 -v' \\<^ = e (36) 

A m in 

and the proof is complete. 

(ii) It will be shown that if ««, i R then there exists a G > such that 

TFM + Ov x = b„ , (37a) 

TF(th) + Gv, = b v , (37b) 
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and 



i'« r* Pi • 



(37c) 



In (37a) v M denotes the terminal vector corresponding to an input 
approaching the constant vector &„. The solution of the network equa- 
tion (3) with b = &„«, and initial condition i>, is v(t) = v,. 

The proof is by simply observing from the definition of R and Lemma 1 
that there exist Vi and G > such that 



T{F(v a ) - F(v')} +G(i; w -v') = 0. 



(38) 



At least when the static diode characteristics are of the conventional 
exponential type the region R is easily obtained as shown below: 



(v' - vYT{F(v') - F(v)\ = v f TD v .F(v) 



(39) 



where 



and 



V = v — V 



D„. = diag (exp X,y( , exp \ 2 v' 2 , ■ ■ ■ , exp X 2n y2„). 
The necessary and sufficient conditions for the positivity of v'TD„-F(v) 




► V 2L-1 



Fig. 4 — Region R for the z'th transistor (exponential diode characteristics). 
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are, from (19), 



, K X 2 ,-_i thi-i exp (X 2l --it>a,--i) < J^ 

tl ~ X 2 , ' a 2 . exp (X 2 ,v 2l ) ~ «2 



for * = 1, 2, • ■ ■ , n. (40) is equivalent to 



X 2 , «2 



111 



a 2 i«i 



^ (X„-it&_, - X2,y 2 ,) ^ In 



(hi 



[flii-li \fl2i-l0l2) J 

for i = 1, 2, • ■ • , n. The region R is plotted in Fig. 4. 



(40) 



(41) 
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